Understanding destruction of nth-order quantum coherence in terms of multi-path 

interference 
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The classic example of the destruction of interference fringes in a "which- way" experiment, caused 
by an environmental interaction, may be viewed as the destruction of first-order coherence as defined 
by Glauber many years ago jl[ . However, the influence of an environment can also destroy the nth- 
order quantum coherence in a quantum system, where this high order coherence is captured. We refer 
to this phenomenon as the nth-order decoherence. In this paper we show that, just as the first-order 
coherence can be understood as the interference of the amplitudes for two distinct paths, the higher 
order coherence may be understood as the interference of multiple amplitudes corresponding to 
multiple paths. To see this, we introduce the concept of nth-order "multi-particle wave amplitude". 
It turns out that the nth-order correlation function can be expressed as the square norm of some 
"multi-particle wave amplitude" for the closed system or as the sum of such square norms for the 
open system. We also examine, as a specific example, how an environment can destroy the second 
order coherence by eliminating the interference between various multiple paths. 

PACS numbers: PACS number(s): 03.65-w, 32.80-t, 42.50-p 
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I. INTRODUCTION 

A most profound concept in quantum physics is quan- 
tum coherence. The first order version of quantum co- 
herence can be directly manifested by the superposition 
of two quantum states. This is quite similar to the opti- 
cal coherence in Young's double experiment . From the 
standpoint of the photon- detection theory by Glauber, 
this simple coherence can be mathematically depicted by 
the first order correlation function However, in quan- 
tum mechanics, this first order coherence phenomenon 
does not sound very marvellous since the same circum- 
stance can also occur in a classical case, such as the op- 
tical interference in the above mentioned Young's dou- 
ble experiment. In fact, by the the first order corre- 
lation function only, it is impossible to distinguish the 
natures of a laser light field and a conventional light 
field with identical spectral properties. As an effective 
remedy, Glauber's n-th order quantum correlation func- 
tion, which accounts for various intrinsic ( n-th order) 
quantum coherence effects, such as the intensity-intensity 
correlation measurement in the Hanbury-Brown-Twiss 
experiment^], was introduced jjj. Indeed, this function 
reflects the intrinsically quantum features of coherence 
beyond the classical analogue. 

The quantum coherence embodies the wave nature in 
the world of microscopic particles. On the other hand, it 
is very fragile and can easily be destroyed by a "which- 
path(way)" experiment ^, Q], or by an environmental 
interaction. This phenomenon of destruction of coher- 
ence is usually referred to as quantum decoherence || |) . 
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The phenomenon of losing the coherence described by 
the first order correlation function is defined as the first 
order quantum decoherence. Motivated by the consid- 
erations in fundamental quantum measurement problem 
0, ||, and also by the attempts to preserve quantum 
coherence of qubits in quantum computing ^ |l(j, pd| , 
many recent experimental and theoretical investigations 
have been focused on revealing the physical mechanism 
of the decoherence problem, e.g., see [|2j. According to 
these studies, this first order decoherence can be roughly 
understood through the quantum entanglement of the 
considered system with the environment or the measur- 
ing apparatus. 

Obviously, this entanglement implies a "which- 
path(way)" detection [|| in the single particle picture. 
Precisely speaking, in an initial coherent superposition 
\ip s ) = ^2c n \n) , each system state |n) corresponds to 
a "path" and many "two path" interferences are re- 
flected in the square norm of the spatial wave function 
(x|'0s)-Thus by considering (x\ip s ) the quantum coher- 
ence can be captured to some extent. After the interac- 
tion, each "path" is correlated with an environment state 
|e ra ) to form an entangling state |i/>t) = c n\ n ) |e n ). 
Here, the different states |e„) distinguish among the 
"paths" of different |n) and thus record the information 
of each "path" . The interference terms in the spatial 
intensity I(x) — Tr((x\?pT} (V'tI^)) will disappear when 
the environment states |e„) are completely distinct, i.e., 
(e m |en) = S mn . In that case each path is labeled by an 
environment state . 

The above well-known explanation of the first order de- 
coherence in terms of "which-path(way)" detection mech- 
anism is simple but very profound. However, it is not yet 
clear whether this mechanism can be used to elucidate 
the nth order quantum decoherence (n-QDC), the de- 
struction of quantum coherence described by Glauber's 



2 



nth order correlation function. The difficulty is we do not 
exactly know what are the "paths" and the corresponding 
"which-path" detection. Most recently, we have touched 
the second order quantum decoherence problem ||l3| by- 
passing this difficulty. As a matter of fact, in our treat- 
ment, we did not define the concept of "path" directly. 
The concrete calculation in the ref. motivated us 

to consider the "which-path" picture of the higher order 
quantum decoherence in general. 

In our present investigation, it is crucial to show, we 
notice that, for a close system, in some cases Glauber's 
nth-order correlation function is the square norm of the 
nth-order "multi-particle wave amplitude" , which will be 
defined later in this paper, while for an open system, it 
can become a sum of the square norms of the nth-order 
"multi-particle wave amplitude" over the states of an en- 
vironment or an apparatus interacting with this open sys- 
tem. This observation is crucial in our present investi- 
gation. As an effective wave function, this multi-particle 
amplitude can be shown to be a supposition of many 
generalized "paths" (the multi- particle paths or simply 
multiple paths). Thus, the higher order coherence may be 
understood as the interference of multiple particle ampli- 
tudes. With this conception a generalized "which-path" 
detection may be established in terms of "multi-particle 
paths" as the physical mechanism of higher-order deco- 
herence for some examples. 

In section 2, we will briefly explain Glauber's nth or- 
der quantum coherence in terms of the sin gle and multi 
photon effective wave functions used in [[14], [15) . As their 
generalizations, in section 3 the concepts of multi-particle 
path and multi-particle wave amplitude are introduced 
for both close and open bosonic systems. Especially, 
the nth order correlation functions of bosonic systems 
will be studied. In section 4, an interacavity model with 



two bosonic modes is used to demonstrate 2-QDC as a 
"which-path" detecting process. In section 5, the exact 
solution obtained in the appendix is utilized to show the 
dynamical process of 2-QDC , which is caused by the en- 
tanglement with the environment or an apparatus and 
indeed can be explained as a generalized "which-path" 
measurement for the explicitly -defined multi-particle 
paths. 



II. THE nTH-ORDER COHERENCE FOR 
QUANTIZED LIGHT FIELD AND 
MULTI-PHOTON WAVE AMPLITUDE 

In quantum mechanics, a pure state is a superposi- 
tion — ^2 k Ck\k) of many components |fc) , but the 
corresponding mixture p — \ck\ 2 \k)(k\ can describe 
the same classical probability distribution \cu\ 2 ■ However, 
and p represents different quantum realities. Usu- 
ally it is said that the two components of a pure state is 
more coherent than those of a mixed one. This coherence 
property is obviously reflected by the intensity interfer- 
ence of two "paths" corresponding to the two components 
in the quantum state. For the system of one particle , 
only single particle property is relevant for this obser- 
vation. In this sense, only the intensity interference ex- 
periment is essential for one particle system. For many 
particle system, however, there exist many experiments 
(such as the Hanbury-Brown-Twiss experiment and the 
intensity-intensity correlation measurement) to show the 
much richer natures of quantum coherence. 

In order to study quantum coherence in many parti- 
cle system, Glauber introduced the so-called n-th order 
quantum correlation function (n-QCF) 
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(1) 



for the electro-magnetic field E a {t) in different mode a. 
Here, E a {t) is the annihilation operator of mode a at 
time t in the Heisenberg picture, E\(f) is the correspond- 
ing conjugate operator, and the density matrix p rep- 
resents the initial state of many-mode electro-magnetic 
field. Of course, this formalism can also be used to study 



the coherence property for any quantum many-body sys- 
tem. 

Furthermore, to describe the higher order coherence, 
Glauber also defined the nth-order coherence function in 



the form |16 



9a/3 = g[ai,a 2 , • • • , a n ; (3i,02, ■ ■ ■ ,Pn) = 
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(2) 



Obviously, this coherence function is defined by a ratio of the off-diagonal elements of the reduced multi-time 
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density matrix g — {g a p : a = (ai, a>2, • • • , = 
{fili &2> ■ ■ ■ > /?n))to its diagonal ones. It is easily observed 
from this definition that these off-diagonal elements rep- 
resent the coherence effect, and each one g a p correlates 
two diagonal ones g aa &nd gpp. Thus, we can understand 
the coherence function as measuring the degree of coher- 
ence for the two diagonal elements, which correspond to 
two different "paths" . From this point of view, there is 
obviously no quantum coherence for a completely mixed 
state with vanishing off-diagonal elements. In the follow- 
ing discussion, we will give simple examples to illustrate 
this viewpoint. 

Now let us turn to an instance |Q from which we 
can see clearly how the higher order quantum coherence 
effects given by 2-QCF is revealed in the multi-particle 
picture. For a single photon, the coherent superposition 
\s) = ^(|lfe) + |l-k}) °f t w0 states with opposite wave 
vectors k and k' possesses the first order quantum coher- 
ence which can be described by the interference fringes 

G w (r,r,t) = (s\E-(r,t)E+(r,t)\s) 

= \(Q\E+(r,t)\s)\ 2 occos 2 (/cr), (3) 

where 

E + (r, t) — E k a k exp(ikr — icOkt) 

is the photon field operator with positive frequency for 
the annihilation operator a k . It should be noticed that 
the diagonal element £?W (r, r, t) of the first order correla- 
tion function is just the square norm of the single photon 
wave-packet |I5l 

{0\E+(r,t)\s) = -±={{0\E+(r,t)\l k ) + {0\E+(r,t)\l k )} 

= —(E k e ikr - iUht + £_ fc e~ lfcMW - fc *)(4) 
V2 

Therefore, (r, r, t) represents the interfer- 

ence between the "two paths" (0\E + (r, t)\lk) and 
(0\E + (r,t)\lk)- This just developes the corresponding 
concept in Young's double experiment for quantized light 
field. The above reformulation of the first order quantum 
coherence implies that "two paths" are necessary for the 
interference phenomenon. 

If we consider the two particle state |lfc, l k >) with a 
single component, an interesting situation arises where 
the first order quantum coherence does not appear, but 
we can see the second order effect through the second 
order quantum correlation function 

G (2) (n, r 2 ,h,t 2 ) = 2E 4 k {l + cos[(fc - k')( n - r 2 )]}. (5) 

Unlike the case of first order coherence, in this case there 
do not appear the obvious two or many "paths". Nev- 
ertheless, the interference phenomenon can still be cap- 
tured in a similar way with the introduction of gener- 
alized "path" . We understand the generalized "path" 
as described by The "two time" correlation function 



G( 2 )(ri, r 2 , ti, tz) = |V>| 2 where ip is the "two-photon wave 
function" 

V; = ^(ri ) r2,ti,t2) = (00|£ + (r 2j t2)^ + (ri,t 1 )|lfc,lfci> 

_ ^2 e -2iw fc t jgifcri+ifc'r'2 _|_ e ifc'ri+ifer 2 j ^ 

if) was invoked as a two photon effective wave-function, 
and it was also called the biphoton wave packet for the 
photon field E + (r, t) . Especially, we remark that 
the biphoton wave packet tp is a coherent superposition of 
two "probability amplitudes" corresponding to two "two- 
photon paths" 

(00|£+(r 2 ,i 2 )|O fc , lfci)<O fc) l fc /|S + (ri s *i)|l fcj 1*) 

and 

(00|^ + (r 2 ,t 2 )|I fc ,O fe O<ffc,0^|S+(r 1 ,i 1 )|I fe ,I fe O 

Starting from Glauber's standpoint and proceed- 
ing along, we come to the conclusion that a set 
{G( n )[ai,£i;a; 2 ,i 2 ; ••• ;a„,t„] \n = 1, 2, ...} of correlation 
functions, rather than a single one, is indispensable to de- 
scribe comprehensively the wave-particle dual nature in 
the quantum world of many particle system. The above 
example of second order coherence shows that the con- 
ception of two-path interference in single photon picture 
still works with a proper generalization of the concept 
of "path". Hence it is quite natural to seek a general- 
ized "which-path(way)" measurement as the mechanism 
of higher-order quantum decoherence. 

III. MULTI-PARTICLE AMPLITUDE FOR 
FREE BOSONS 

In this section the conception of two photon effective 
wave-function will be generalized. It will be applied to 
the study of general quantum systems of identical parti- 
cles. We first discuss the spatially-homogeneous case for 
the sake of simplicity. 

We consider a homogeneous bosonic field with two 
modes \V) and \H). The generalized field operator in 
time-domain 

4>(t) = c v b v e- iuvt + c ff b H e- iuJHt 

= c v (t)b v + c H (t)b H (7) 

is an annihilation operator with respect to the superpo- 
sition state 

\+)=c* v \V)+c* H \H). (8) 

Here , bn and by are the annihilation operators of the 
boson system ; cy and ch satisfy the normalization re- 
lation \cy\ 2 + \ch\ 2 — 1. Without loss of generality, we 
take cy = ch — l/v2 . This means we consider the 
measurement to detect the polarized boson along the 45° 
direction in the V — H plane. We call <j> a "measuring" 
operator . 
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Corresponding to 0, the generalized second order cor- 
relation function (ll|] 

G (2) - {IvIh\^ (h)^ (t 2 )4>(t 2 )^{ti)\lv Ah) 

= KO,0|^(t 2 )^(ti)|lv,l ff )| 2 = |*(ti,ta)| a . (9) 

The two time wave function 

*(ti,*a) = (O,O|0(t 2 )0(ii)|l y ,l ff ) 

can be understood in terms of the two "paths" picture 
from the initial state |ly, 1h) to the final state |0,0): 



h i \ Cff(t P 


|iv,o ff ) c ^ 2) 10,0) 


\ cy(tl) 


Cff(t 2 ) y 




|0y,l ff ) 



The two "paths" are just associated with the two ampli- 
tudes forming a coherent superposition 

*(ii,*2) = eye// exp(— zwyi2 — iVff) + 

c/fcy exp(— iu H t 2 — iu v ti) (10) 

Correspondingly, the second order correlation function 

G< 2 > = 2\c v c H \ 2 {l + cos[(lu v - uj H )(t 2 - h)} (11) 

The above discussion for the second order quantum 
coherence is applicable to the higher order case. Our 
arguments in this paper are based on two novel obser- 
vations: a. The generalized field operator = ^2c n b n is 



specified for a quantum measurement about a superposi- 
tion single particle state \(f>) = c^\n) . b. For a certain 
initial single component state \sq) of N particles system, 
the n-th order quantum correlation function 

G (B) (ri,rs,,---,r n ,ti,t2,---,*n) = lV' (n) | a (12) 

can be written as the norm square of an effective wave 
function ip^ n \ which is just a superposition of many am- 
plitudes. 

Let us consider the third order situation as an example. 
Let the initial state be |lfj,2y). Then the generalized 
third order correlation function 

G<® = {2 v lH\^Hti)^Ht 2 )^\h)^{t^{t 2 )^{t l )\2 v ,l H ) 

= mo\4>(t 3 )j>(h)kti)\2vAH)\ 2 

= \y(h,t 2 ,t 3 )\ 2 (13) 
is a norm square of the two time wave function : 

*(ti,t 2 ,t 3 ) = {o,o\$(t 3 )4>(t 2 Mh)\2 v ,i H ) 

= V2c v 2 c H e- Wv{t3+t2) ~ lUHtl + 
V2c H c v 2 e' WHt2 ^ v{t3+tl) + 
V2c H c v 2 e- WHt3 - Wv{t2+tl) (14) 

Each term in the above effective wave function is con- 
tributed by the corresponding one of the four "paths" 
from \2vAh) to |0,0): 



\2 v ,0h) 


Cv ^ \1v,0b) 








cv(*3)\ 


\2v,1h) ^ ) \1vAh) 


C ^ 2) llv.Oir) 


c ^l ) |0,0) 




\Cv(t2) 


c H (t 3 ) y 




|0y,lff) 





In terms of the effective 3-time-wave function defined 
above, the third order correlation function is explicitly 
written down: 

G (3) = 4\c v 2 c H \ 2 (^+cos[(uj v -Lj H )(t 2 -ii)] + 

cos[(lu v ~ v H )(t3 - h)] + cos[(wy - LU H )(t 2 - t 3 )]) 

It shows the quantum interference in the time-domain. 

The above analysis is valid only for the case where 
the considered system is isolated from an environment 
and not measured by a detecting apparatus-a detec- 
tor. For our purpose, we need to consider an open 
system S interacting with an environment (reservoir) 
or a detector E, and we must extend the concepts of 
multi-particle(time)-wave functions and the correspond- 
ing many-particle paths defined above. To do the gener- 



alization , we first invoke the effective field operator 

B{t) = U\t,0)4>(0)U(t,0) 

= exp(iVt)$(t) exp(-iVt) (15) 

where 0(0) = cyby + cnb H has been given by Eq.(||). 
Then instead of the free time evolution governed by 
the free Hamiltonian Hq, we use the evolution operator 
U(t, 0) governed by the total Hamiltonian 

H = H + W + H E = H + V 

taking into account the role of the interaction between E 
and S . Here, He is the free Hamiltonian for E. If we only 
consider an ideal quantum decoherence process without 
dissipation, V should possess the nature of quantum non- 
demolition: [H , W] = |@, H. 
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Let the states |n) = |ny,n^) be the common eigen- 
states of H and W corresponding to the egen-vaslues 
E n and V(n)(n — (ny ,nn))- If the initial state of the 
total system is 

|V(0)) = \d> a ) ® (16) 

where |0 S ) and are some specially-given initial states 
of S and E respectively, we can define the effective two- 
time state vector 

\iMt,t'))=B(t')B(t)\m)- (17) 

as a reasonable generalization of the effective "two-time 
wave function" given above. Its norm is just the second 

I 



Here, the summation ranges over the complete set of 
states 1/3) of E, and each term in the sum is a norm 
square of the effective two particle wave function 

*/j(*i,ta) = (0,P\B(t')B(t))\l H ,lv,<f>E) (20) 

for the open system. 

From the above calculations for the second and third 
order quantum decoherence, we observe that for a 
specially-given initial state, a higher order correlation 
function may be explicitly written down as the norm 
square (or its sum ) of the multi-time-wave function, 
which is a coherent superposition of several complex com- 
ponents associated with the generalized many-particle 
paths. It is pointed out that this kind of many-particle 
path is not a simple-product of single-particle paths, but 
it can be determined by the specially designed measure- 
ment. 



IV. GENERALIZED WHICH-PATH 
DETECTION IN AN INTRACAVITY MODEL 

In this section, an intracavity model is presented to 
demonstrate in multi-particle picture the "which-path" 
detection associated with higher-order quantum decoher- 
ence . 

In a recent paper E3], we have studied the problem 
of 2-QDC for a cavity-QED system. The concrete cal- 
culation in the rcf. (Tj| shows that the 2-QDC effects can 
indeed be observed in the proposed experiment. But it 
involves dissipation effect losing energy. However, it is 
well known that quantum decoherence can still occur for 



order correlation function 

(<Mi,OIV^(M')} = Tr{mB\t)B\t')B{t')B(t)) 
= G^[t,t',p(0)}, (18) 

for the density matrix ,5(0) = |<0(O))(^(O)|. 

In fact, due to the non- demolition interaction not re- 
sulting in dissipation, the basic dynamic properties of the 
open system do not change even in the presence of E . 
If we choose \4> s ) = \1h, ly), there are only two "paths" 
from the initial state |ly, 1h) to the final state |0,0) for 
S , and 1 0,0) is the unique state which can be reached 
by the action of B(t')B(t). Then, 



(19) 



an energy conserving system. So in principle dissipa- 
tion is not indispcnsible for the discussion about deco- 
herence effect. For this reason it is natural and inter- 
esting to consider pure decoherence process without dis- 
sipation. The pure decoherence can be well understood 
through the quantum entanglement of the considered sys- 
tem with the environment or the measuring apparatus. 
For a model with pure decoherence process, many con- 
cepts (such as multi-particle path and the correspond- 
ing which-path detection ) can be made much clearer. 
Unlike the approximately-solvable model treated in the 
ref . Jl3| , which loses its energy and coherence simultane- 
ously, the model proposed in this section, a bosonic sys- 
tem of two modes interacting with an external system of 
many harmonic oscillators, does not dissipate its energy. 
This property makes the model exactly solvable, and as a 
result the problem of higher-order quantum decoherence 
can be studied in a straightforward way. 

By taking h — 1, the model Hamiltonian H = Ho + V 
is defined by 

H = uyblbv, (21) 
^ = + S }2,[dv{u 3 )b\,b v +d H {ujj)tf H bH\ 

3 3 

a]+aj), (22) 

where Ho is the free Hamiltonian of the system, V the 
free Hamiltonian ■ LUjOjCij of the reservoir (or a detec- 
tor) plus a non-demolition interaction between the sys- 
tem and the reservoir; and Sy(6y), V h (J)h) the creation 
(annihilation) operators for the two modes with frequen- 
cies ujy 0) and ojh = 0. The operators a^{a,j) are 



{i> B (t,t')\Mt,t')) = X>B(M0k/3>MV>B(M')> = ]T|<o,/3|?MM'))l 

n,/3 

= J2 \(0,(3\B(t')B(t))\l H , Iv^e)\ 2 - £ |tf/9(ti,*2)| a 

P 
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creation (annihilation) operators of the reservoir modes 
of the frequencies u>j . The frequency-dependent constant 
dnitOj) (dv{<jJj)) measures the coupling constant between 
H (V) mode and the j mode of the reservoir. 

The most important feature of the model is the non- 
demolition condition [Ho, V] = 0. It means that the sys- 
tem does not dissipate energy to the reservoir. On the 
other hand, the system can leave imprint on the reser- 
voir since, for different number states \ny, ur), there are 
different interactions 

y~][n v d v (u)j) + n H d H {uj)]{a] + dj) 

3 

acting on the oscillator reservoir with different driving 
forces ~ nydy(u)j) + n#dff (wj). When there is only one 
mode in the external system(reservoir) the whole system 
can physically be described by an intracavity model: Two 
mode cavity field interact with a moving wall of the cav- 
ity, which is attached to a spring and can be regarded 
as a harmonic oscillator with a small mass [4,8]. The 
fields are coupled to the cavity wall (a moving mirror) by 
the radiation pressure forces in proportion to the photon 
numbers V H bu and b v by. 

For the above introduced model, we can discuss the 
higher order decoherence problem in the Heisenberg pic- 
ture by explicitly defining the many-particle "which - 
path" measurement. The second order coherence is di- 
rectly determined by the second order correlation func- 
tion. 

G[t,t',p(0)} =Tr(p(0)BHt)BHt')B(t')B(t)) (23) 

which is defined as a functional of the density operator 
p(0) of the whole system at a given time 0. Here, the 
bosonic field (measuring) operator 

B(t) = cxp(iVt)(cHbH + cyby exp(-iujyt)) cxp(-iVt) 

(24) 



is defined for the interacting system. Like the operator 
defined for the non-interacting system, it also describes 
a specific destructive quantum measurement [4] with re- 
spect to the polarized states 

\+)=c* H \H)+c v \V) 

where ch and cy satisfy the normalization relation 
\ch 1 2 + \ c v\ 2 = 1- Without loss of the generality, we take 
ch = cy = l/\/2, considering a specific measuremeant. 

To examine whether the macroscopic feature of the 
reservoir causes the second order decoherence or not , we 
consider the whole system in an initial state 

|V(0)) = |WW®l{0j}>, (25) 

where |{0j}) is the vacuum state of the reservoir. Here, 
we have denoted the general Fock states of the many 
mode field by |{rij}) = \ni,n 2 , ...). Because the present 
discussion concerns the external system interacting the 
considered system , the conceptions presented in last sec- 
tion must be alternated. Actually, in stead of the effective 
"two-time wave function" , we use the effective two-time 
state vector 

\^ B {t,t'))=B{t')B{t)\^)). (26) 

Then we can re-write the second order correlation func- 
tion as 

G[t,t',p(0)} = (^B(t,t')\^ B (t,t')) (27) 

It is interesting that the effective state vector can be 
evaluated as the superposition 



l^ B ( M ')) = I e *mo )t ' [expH ^ 

, (28) 

of two components for the two paths from the initial \0h,0v) ■ It should be noticed that the effective actions 
two particle state \1hAv) to the two particle vacuum of the reservoir 



V(m, n) = ^ Vj(m, n) = ^ ujja}-aj + '^^(dy(uij)m + da (wj)n)(aj + dj) (29) 



can label the different paths and record the path informa- tion in the reservior. Thus , this generalized "which-path 
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measurement " leads to the second order quantum deco- 
herence. 

The above result clearly demonstrates that, with the 
presence of the reservoir, the different probability ampli- 
tudes (~ exp(— iuiyt') and exp(— iuiyt)) from ly) to 
|0j?,0y) entangle with the different states 

(i e iv(o,o)t' e -iy(i,o)t' e iy(i,o)t e -iy(i,i)t|| _,^ 

and 

i e i^(o,o)f / e -»V'(o,i)t' e iv(o,i)t e -it>-(i,i)t^o i })) 

of the reservoir. This is just the physical cause of the 
second order quantum decoherence. In the following sec- 
tion an explicit calculation of the second order correlation 



function will be given to illustrate this crucial observa- 
tion. 



V. DYNAMIC DECOHERENCE IN HIGHER 
ORDER CASE 



After a straightforward calculation the second order 
correlation function can be expressed in a factorization 
form [9]: 

G[t,t',p(0)} = -[l + Me^"*') n*»] ( 30 ) 



where each factor 



J 



Fj = (0j\e 



^•(iA)t P -i^(o,i)t„iVi(o,i)i' -iV5(i,o)t' ^(l.o^p-i^a.iHin.s = 



\0j) = (0,1^(^)10,) 



(31) 



is a two-time transition amplitude of the j'th mode of 
the reservoir. Obviously, the term FT . Fj measures the 
extent of coherence and decoherence in the second order 
case. It plays the same role as the decoherence factor of 
the first order decoherence [Q . So it is also called the 
decoherence factor. 

In the following, to give the factor Fj explic- 
itly, we adopt the Wei-Norman method 0, |1| 
to calculate the effective time evolution defined by 
Uj(te). It can be imagined as an evolution governed 
by a discrete time-dependent Hamiltonian H(t) dom- 
inated by Vj(l, 1), -Vj{\, 0), Vj(l, 0), -^(0, 1), fyo, 1) 
and — Vj(l, 1) in six time-intervals [to — 0, ii = 
t\,[t u t 2 = 2t],[t 2 ,*3 = 2t + t'],[t 3 ,U = 2t + 2t'],[U,t 5 = 
3t + 2t'}, [t$, to = 4t + 2t'\ respectively. In the fc-th step 



of calculation, we take the final state of (k— l)-th step as 
its initial state. Therefore, we obtain iij(te) as the sixth 
step evolution 



(f 6 ) = e 9?3( t 6) a i + e 92j( t 6)a 3 t aj e S3j( t 6)ai e sl3(*e) (32) 



Here, g%j{t^){k = 1,2,3,4) are the coefficients that can 
be explicitly obtained. But for the calculation of the j'- 
th component 



Fj — exp[<7|(t 6 )] 



(33) 



of the decoherence factor, we only need to know g\y The 
detailed discussion in the appendix gives 



J 



9%{h) = -\[d H { 3 )-d v {j)f^[\uj (t'-t)] + ^[dUj)~dl{j)] 

[ojj(t f — t) + 2(1 — cos(wj[t' — t]) smuijt] + (1 — 2 cosojjt) sin(u)j[t' - 
= -Rj(t -t') + «%(t,t') 



r 



(34) 



It is noticed that the real part 

-Rj (t - t') = - A [ dH (j) -dv(j)} 2 sin 2 [L, (f 



t)] 



only depends on the time interval t'—t, but the imaginary 
part £lj(t,t') depends on both t' and t' as a two time 
function. 



Because —Rj can not exceed zero, the norm \Fj\ = 
e~ R iWoi the factor Fj can not exceed one . Then from 
the arguments about the first order decoherence in our 
previous works on quantum measurement theory [frl , it is 
concluded qualitatively that the factorization structure 
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of the decoherence factor 

N N N 

F = Y[ Fj = H e^C*.* 7 ) • H e -^(*-*') 

2=1 2=1 2=1 

= \F\exp[ifl(t,t')} 

implies the vanishment of the second order correlation 
in the macroscopic limit that the number N of particles 
making up the reservoir approaches the infinity. This is 
because 

N 

\F\=exp[-Y,fyM]^0 

2=1 

as N — > oo since Y^j=i Rj ^ s a diverging series or 
a monotonously-increasing function of t — t' under some 
reasonable conditions. 

In order to demonstrate the above conclusion quan- 
titatively, we give the numerical results for the second 
order decoherence for different numbers N of the quan- 
tum oscillators. As N increase, these results are illus- 
trated in Fig.l. In the numerical calculation, the coupling 
constants {dy(j)} take random values in the domain 
[0.8, 1.0], the coupling constants {du(j)} in [0.2, 0.4], and 
the frequencies {ujj} in [0.5,1.5]. The other parameters 
are given in the caption of the figure. 

From both the above graphic illustration given by the 
numerical calculation and the analytic result given by 
Eqs.d^HHJ) , we see clearly that the second order cor- 
relation function depends not only on the time interval 
t' — t, but also on t or t' , In fact it follows from the 
analytic result 

Imie^-^YlFi) = ^\F\cos[uj v (t-t')+Sl(t,t')] 

(35) 

that the amplitude of the second order correlation func- 
tion is mainly determined by the time interval t' — t, 
but the phase is determined by the two time parame- 
ters . What is more important is that as the number of 
the quantum oscillators increases, the second order co- 
herence vanishes faster and faster, and the amplitude for 
quantum revival becomes smaller and smaller. With a 
reasonable extrapolation, it can be predicted that, when 
the number of the quantum oscillators approaches the 
infinity in the macroscopic limit, the second order coher- 
ence will vanish in a very short time and thus no quantum 
revival phenomenon can be observed. 

VI. CONCLUDING REMARKS 

In this paper we first depict the high order quantum 
coherence of a boson system by introducing the concept 
of multi-particle wave amplitude. For some cases with 
the specifically- given initial state , we show that the 




0.2- 

10 20 30 40 t 

g d 

lr 




ltT ' ' 20 ' ' 30 ' ' "~40 t 



FIG. 1: The horizontal axe denotes time period t' — t, 
the vertical axe denotes the second order correlation func- 
tion G = G[t, t' , p(0)], parameters ujy = 1.0, (a)iV = 5,t — 0, 
(b)N = 5, t = 5, {c)N = 10, t = 0, (d)N = 10, t = 5. 
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norm square of multi-particle wave amplitude (or a sum 
of the norm squares for an open system) gives the high or- 
der correlation function. As an effective multi-time wave 
function, this amplitude can be shown to be a superpo- 
sition of several "multi-particle paths" . When the envi- 
ronment or an apparatus entangles with them to form 
a generalized "which-path(way)" measurement, the high 
order quantum decoherence happens dynamically. Some 
explicit and general illustrations are presented in this pa- 
per to construe our observation. But to prove our conjec- 
ture that any high order correlation function indicating 
the existence of high order quantum coherence can be ex- 
pressed as the norm square of a properly-defined effective 
wave amplitude for many particle system, there is still a 
long way to go. If this conjecture is true under certain 
general conditions, then there still exist the problem of 
clarifying these conditions. Moreover, experimental pro- 
posals that can be implemented at least in principle are 
still unavailable. 

Our present investigations shed a new light on the un- 
derstanding of quantum coherence . According to argu- 
ments in this paper the high order quantum coherence 
can also be sculptured as the generalized interference 
phenomenon by two "multi-particle paths " , and the in- 
trinsically quantum features (higher order ones) of coher- 
ence beyond the classical analogue reflected by the spatial 
interference of two paths in classical electromagnetic field 
can be theoretically unified in a framework to embody 
the wave - particle duality in the quantum world. There- 
fore, as the essential element, the "which-path(way)" de- 
tection in both the original and the extended versions, 
naturally provides a complete decoherence mechanism in 
understanding quantum measurement and the transition 
from quantum to classical mechanics. 
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tk ~ ife-i take different values in the six different steps: 



A 1 


— R 1 — 

— u 'j ) Pj — 


■y 1 - 


ay [ujj ) -t au [uij ) , J. i — 


,x' 2 

a / 


, , «2 

- ~Uj,Pj 


- ~? 


— av\y>j), J-2 — t, 


a? 


= ">3>0j = 




d v {uj),T 3 = t', 




= 


= 7* 


= -d H (oJj),T4 = t', 




= o>j>$ = 


7? = 


d H (ujj),T 5 = t, 






% 6 = 


-d v (u)j) - d H (uj),T 6 



Because of the fact that the four operators rij = a^dj, 

Oj , dj and 1 form a closed algebra - the Heisenberg-Wely 
algebra, the unitary time evolution operator at each step 
takes the following form (Wei-Norman theorem) 

fi*(T) = e^ & U^ T ^e&™e^ T \ (38) 

for T £ [tk-i,tk] in a special sequence. Here coefficients 
9sj(T){s — 1, 2, 3, 4) are functions of T to be determined. 
The benefit of the above form is that only the coefficient 
g%j (T) is needed in the calculation of the average value at 
the vacuum state. So we can largely reduce the complex- 
ity of our calculation as we need to pay attention only to 
things concerning g\^(T) . 



Substituting (T) into the Schrddinger equation 
d 



i—vlt = h)u k p 



dT 3 



we find the coefficients g^j(T)(s = 1,2,3,4) satisfy the 
following system of equations: 



d k k 

Using the results 
d 





d k 

dT 92j 


= 


-9ij 


d k 
dT 9 2j 


= -i® 




d k 
dT 93j 


= ~H 




d k 
dT 93j 


= 



(39) 



APPENDIXES 



-^g\ d = -ia)g\ 3 - i# , —g% = -i^g^ (40) 



In this appendix, the Wei-Norman method JT^, |l8| is 
adopted to calculate the second order decoherence factor 
Fj. The calculation is completed in six steps. 

During the time period [tk-i,tk](k — 1,2, •••,6), Let 
Wj(t,t') be an time evolution dominated by the single 
particle Hamiltonian 



hj = aj&jdj + /^aj + rfcLj 



{k = 1,2, ■••,6}. (36) 



The coefficients oft,0j, 7^ and the time intervals Tk 



obtained by simplifying the above system of equations , 
we get the solution 

9ij(T) = (g^tk-x) + ^e-^^-O _ 

g%{T) = gi j (t k _ l ) + \{g^{t k . l ) + %) 



(e -i«**(T-t fc _0 + 



{41) 
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Notice that, to obtain the above result we have used the h \ fc_i/, . li r fc-i /, . ft 

initial conditions 



<&(*k-l) = ^W, (42) (e -i-jT t _ 1)+i W < (44) 

^(tfc-i) - ^-'(tfc-i) (43) 

for each step and the initial conditions <7y(to) = 
ffSjC^o) = for the first step. Then, we obtain a set 

of iteration equations T , . . ,. ... ...... 

Iterating six times with ditterent initial conditions 

ok ok and coefficients, the final result of g\{to) is obtained as 

J J a . on 



[1] R.J. Glauber, Phys.Rev.130, 2529(1963); 131, 

2766(1963). [11 
[2] H. Hanburg-Brow and R.Q. Twiss, 

Phil.Mag.45,663(1954); Naturel78, 1046(1956); [12" 

Proc.Roy.Soc.A242,300(1957). 
[3] S. Durt, T. Nonn, and G. Rampe, Nature 395, 33(1998). 
[4] E. Buks, R. Schuster, M. Heiblum, D. Mahalu, and V. 

Umansky, Nature 391, 871(1998). 
[5] W.H. Zurek, Phys.Today, 44(10), 36 (1991). [13; 
[6] M. Brune, E. Hagley, J. Dreyer, X. Maitre, A. Maali, C. 

Wunderlich, J. M. Raimond, and S. Haroche. Phys. Rev. [14 

Lett. 77, 4887(1996). 
[7] C. P. Sun, Phys. Rev. A 48, 878(1993). C. P. Sun, Chin. [15" 

J. Phys. 32, 7(1994). C. P. Sun, X. X. Yi, and X. J. Liu, [16" 

Fortschr. Phys. 43, 585. 
[8] S. Bose, K. Jacobs, P. Knight, Phys.Rev. A 82, 

3204(1999). 
[9] W.Unruh,Phys.Rev. A 51, 992(1995). 
[10] G. Palma, K. Suominen, A. Ekert, Proc. Roy. Soc. Lon- 



don Ser. A 452, 567 (1996) 

C. P. Sun, H. Zhan, and X. F. Liu, Phys. Rev. A 58, 
1810(1998). 

C. P. Sun,X. F. Liu, D. L. Zhou, and S. X. Yu, Phys.Rev. 
A 63, 2111(2001), C. P. Sun, D. L. Zhou, S. X. Yu, and 
X. F. Liu, Eur. J. Phys. D 13, 145(2001);C. P. Sun,X. 
F. Liu, D. L. Zhou, and S. X. Yu, Eur. J. Phys. D 17, 
89(2001) 

D. L. Zhou, and CP. Sun, |quant-ph/0104038| , LNAL 
preprint (2001). 

M.O. Scully, and M.S. Zubairy, Quantum Optics, Cam- 
bridge University press, 1997, pp 97-129. 
Y. Shih, Adv.At.Mol.Opt.Phys. Vol 41(1999), pp 1-42. 
This definition, which is defined by Klauder and Sudar- 
shan in Fundamentals of Quantum Optics, W.A. Ben- 
jamin, NewYork, 1968, more convenient to describe co- 
herence as interference, has a little difference in the de- 
nominator with Glauber's one. Glauber's original defini- 
tion is 



g[ai,a 2 , ■ ■ ■ ,a n ;/3i,/fc, • • • ,&>] = 



TrlpB^ (ti)B^ 2 (t 2 ) • • • Bl n (t n )B a „ (t n ) ■ ■ ■ B a2 (t 2 )B ai (ti)] 
y/G^MG^M ■ ■■G[a n ,t n ]y/G[/3 1 ,t 1 ]G\j32,t 2 ] ■ ■■G[/3 n ,t n 
I 



[17] J. Wei and E. Norman, J. Math.Phys. 4A, 575(1963). 359(1990). 
[18] CP. Sun and Q. Xiao, Commun.Theor.Phys. 16, 



